Is nonhelical hydromagnetic turbulence peaked at small scales? 
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ABSTRACT 

Nonhelical hydromagnetic turbulence without an imposed magnetic field is considered in 
the case where the magnetic Prandtl number is unity. The magnetic field is entirely due to 
dynamo action. The magnetic energy spectrum peaks at a wavenumber of about 5 times the 
minimum wavenumber in the domain, and not at the resistive scale, as has previously been 
argued. Throughout the inertial range the spectral magnetic energy exceeds the kinetic energy 
by a factor of about 2.5, and both spectra are approximately parallel. At first glance, the total 
energy spectrum seems to be close to fc~ 3 / 2 , but there is a strong bottleneck effect and it is 
suggested that the asymptotic spectrum is fc~ 5 / 3 . This is supported by the value of the second 
order structure function exponent that is found to be £2 = 0.70, suggesting a fc -170 spectrum. 

Subject headings: ISM: kinematics and dynamics — magnetic fields — MHD — turbulence 



1. Introduction 

It is generally accepted that in hydromag- 
netic turbulence the magnetic field tends to 
be more intermittent than the velocity field. 
This is evidenced by many numerical simu- 
lations where the magnetic field is dynamo- 
generated (Meneguzzi, Frisch, & Pouquet 1981; 
Kida, Yanase, & Mizushima 1991; Branden- 
burg et al. 1996; Kulsrud et al. 1997). Further- 
more, linear theory predicts that the exponentially 
growing magnetic energy spectrum increases with 
wavenumber like fc +3 / 2 Kazantsev (1968). Taken 
at face value, linear theory would suggest that the 
magnetic energy spectrum should be peaked at 
the 'resistive' cutoff scale. In the case of the in- 
terstellar medium the resistive cutoff scale would 
be tiny (< 10 10 cm) in comparison with other 
relevant scales (> 10 18 cm). There is of course 
no doubt that there are magnetic fluctuations in 
the interstellar medium at a scale below 10 10 cm, 
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as evidenced by interstellar scintillation measure- 
ments (Goldreich & Sridhar 1995), but it remains 
implausible that magnetic fields at such small scale 
contribute significantly to the energy budget. 

The question of small scale magnetic fields has 
worried theorists for the last decade. Kulsrud & 
Anderson (1992) have shown that, independent 
of the possible presence of magnetic helicity, the 
kinematic magnetic energy spectrum follows the 
Kazantsev /c+ 3 / 2 law, and they speculate that this 
may suppress large scale dynamo action (of a 2 or 
afi type); see also Vainshtein & Cattaneo (1992). 
The case of helical dynamos is now reasonably well 
understood in the case of large magnetic Reynolds 
number. The saturation level of helical dynamos is 
not suppressed, but the saturation time is the re- 
sistive time scale (Brandenburg 2001). This is not 
the result of small scale magnetic fields in general, 
but due to the helical small scale fields that are 
produced by the a effect as a by-product (Field & 
Blackman 2002; Blackman & Brandenburg 2002). 

When the degree of kinetic helicity of the flow 
falls below a certain threshold, no large scale dy- 
namo action is possible and the magnetic energy 
spectrum is peaked at a scale much smaller than 
the forcing scale (Maron & Blackman 2002). This 
may not be so much a concern for stellar dynamos 
where differential rotation is important, but it 
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could be a problem for dynamo action in the in- 
terstellar medium and in clusters of galaxies. An- 
alytic approaches to the nonlinear saturation of 
nonhclical dynamos have corroborated the notion 
that the magnetic energy spectrum peaks at the 
resistive scale (Schekochihin et al. 2002a, b). Nu- 
merical simulations have confirmed a peak well 
below the forcing scale (Cho & Vishniac 2000a; 
Maron & Cowley 2001), but a resolution of up to 
256 3 collocation points is still insufficient to estab- 
lish whether the location of the peak is different 
from the resistive scale. 

So far there has been no evidence for fc~ 5//3 in- 
ertial range scaling of the magnetic energy as in 
the case of an imposed field (Goldreich & Sridhar 
1995; Cho & Vishniac 2000b; Maron & Goldre- 
ich 2001) or as in decaying hydromagnctic turbu- 
lence (Biskamp & Miiller 2000). A difficulty in 
establishing power law scaling is the lack of a suf- 
ficiently long inertial range. In some cases the use 
of hyper-resistivity (where the V 2 diffusion op- 
erator is replaced for example by —V 4 , in order 
to shorten the diffusive subrange) has led to the 
concern that it may cause an artificially enhanced 
'bottleneck effect' with a shallower k" 1 spectrum 
just before the dissipative subrange (Biskamp & 
Miiller 2000). Yet another problem is the possibil- 
ity of a physical bottleneck effect (Falkovich 1994) 
that will be more extreme in three-dimensional 
spectra than in the one-dimensional spectra ac- 
cessible from turbulence experiments (Biskamp & 
Miiller 2000; Dobler et al. 2003). 

In this Letter we use simulations at resolutions 
of up to 1024 3 mcshpoints to study the form of 
the magnetic and kinetic energy spectra in the in- 
ertial range in the case where the turbulence is 
driven at large scales and the magnetic field is self- 
generated. No hyperviscosity or hyper-resistivity 
is used. 

2. Equations 

Here we consider subsonic turbulence in an 
isothermal electrically conducting gas with con- 
stant sound speed c s in a periodic box of size 
2-7T x 2n x 2ir. The governing equations are 
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magnetic field, /j, the vacuum permeability, 

F visc = i/(V 2 u+|VV'« + 2S' VLnp) (2) 

is the viscous force where v — const is the kine- 
matic viscosity, Sy = \{ui^ + Uj^) — hSijitk.k is 
the traceless rate of strain tensor, and / is a ran- 
dom forcing function (see below). The continuity 
equation is written in terms of logarithmic density, 



Dlnp 
Dt 



= -V • it, 
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and the induction equation is solved in terms of 
the magnetic vector potential A, where B = Vx 
A, and 

dA ~ , 

— = u x B + r]V 2 A, 
at 



(4) 



and r] = const is the magnetic diffusivity (or resis- 
tivity); we choose r\ = u, i.e. our magnetic Prandtl 
number is unity. We adopt a forcing function / of 
the form 

f(x, t) = Re{A/ fc(t) exp[ifc(i) • x + i<f>(t)]}, (5) 

where x = (x,y,z) is the position vector, and 
— 7r < <f>(t) < tt is a (J-correlated) random phase. 
The normalization factor is N = foc s (kc s / 'St) 1 / 2 , 
with /o a nondimensional forcing amplitude, k = 
\k\, and St the length of the timestep; we chose 
/o = 0.02 so that the maximum Mach num- 
ber stays below about 0.5 (the rms Mach num- 
ber is close to 0.12 in all runs.) describes 
nonhelical transversal waves with \fk\ 2 = 1 and 
fk = (k x e) / yjk 2 — (k ■ e) 2 , where e is an ar- 
bitrary unit vector. At each timestep we select 
randomly one of 20 possible wavevectors in the 
range 1 < |fe| < 2 around the forcing wavenum- 
ber, k{ — 1.5. 

The equations are solved using the same 
method as in Brandenburg (2001), but here we 
employ a new cache and memory efficient code 4 
using MPI (Message Passing Interface) library 
calls for communication between processors. 

3. Results 

In Fig. 1 we plot magnetic, kinetic, and total en- 
ergy spectra, E M (k), E K (k), and E T — E M + E K , 



where D/Dt = d/dt + u- V is the advective deriva- 
tive, J = V x .B/yUo the current density, B the 



4 We use the Pencil Code which is a cache efficient grid 
based high order code (sixth order in space and third or- 
der in time) for solving the compressible MHD equations; 
http://www.nordita.dk/data/brandenb/pencil-code. 
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respectively, for our largest resolution run with 
1024 3 meshpoints. 5 The magnetic energy displays 
a nearly flat spectrum in the range 1 < k < 5, 
peaks at fc ss 5, and begins to show an iner- 
tial range in 8 < k < 25, followed by a dissi- 
pative subrange over one decade. In the inertial 
range Ey[(k) / Exik) is about 2.5. A plot of mag- 
netic energy per logarithmic wavenumber interval, 
kEM{k), shows that most of the magnetic energy 
comes from a band of wavenumbers around k = 9 
- independent of the Reynolds number once the 
latter is above ~ 400; see Fig. 2. 

The energy ratio in the inertial range is sim- 
ilar to the ratio of the dissipation rates of mag- 
netic and kinetic energies, em/ek, which is about 
2.3; see Fig. 3, which shows that magnetic and ki- 
netic energy dissipation rates approach 70% and 
30%, respectively, of the total dissipation rate, 
£t = £m + £k • The convergence of relative dissi- 
pation rates is compatible with Kolmogorov's con- 
cept of a constant, scale independent energy flux 
across the spectrum, which seems to apply also 
separately for velocity and magnetic fields. This 
picture would be difficult to reconcile if the mag- 
netic energy were to always peak at the resistive 
scale. We emphasize that the inertial range is not 
representative of the total energy which, in turn, 
is dominated by small wavenumbers. The ratio of 
total magnetic to kinetic energies is only 0.4 and 
seems again to be asymptotically independent of 
Reynolds number (Fig. 3). 

By comparing runs at different resolution one 
can clearly see that in the range 3 < k < 20 
the total energy spectrum is shallower than fc~ 5 / 3 
(Fig. 4). This could perhaps be due to the bot- 
tleneck effect that is known to exist also in wind 
tunnel turbulence, where it has been described by 
a weak fc" 1 contribution (She & Jackson 1993). 
There are several reasons why such a bottleneck 
effect might occur. First, recent studies by Doblcr 
et al. (2003) have shown that the bottleneck ef- 
fect is much stronger in shell integrated three- 
dimensional spectra compared to just longitudinal 
or transversal one-dimensional spectra available in 
wind tunnel turbulence. Second, the bottleneck 
effect may simply be stronger for hydromagnetic 



5 These spectra arc, as usual, integrated over shells in k space 
and normalized such that J E^dk = \ {n?) and f EM<ik = 
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Fig. 1. — Magnetic, kinetic and total energy spec- 
tra. 1024 3 meshpoints. The Reynolds number is 

Hrms/Kf) ~ 960. 
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Fig. 2. — Spectral magnetic energy per logarith- 
mic wavenumber interval for runs with different 
Reynolds numbers. The arrow correspond to the 
peak (k = 9) of the four largest runs. 



turbulence due to the dynamo effect that is ex- 
pected to produce magnetic energy preferentially 
at k ~ fcfi?y 2 ;=y 5. ..10, where R c 30 is the 
critical magnetic Reynolds number for dynamo ac- 
tion (Subramanian 1999). Third, numerical effects 
such as hyperdiffusion or other effects causing de- 
partures from the physical V 2 diffusion operator 
could cause an artificial bottleneck effect (Biskamp 
& Miiller 2000; Biskamp et al. 1998), which may 
also explain the extended fc _1 range in compress- 
ible nonmagnetic simulations using the piecewise 
parabolic method (Porter et al. 1998). 

We believe that numerical effects do not cause 
such artifacts in our simulations, because we use 
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Fig. 3. — Approximate convergence of relative 
magnetic and kinetic energy dissipation rates (up- 
per panel) and the relative energies (lower panel) 
as a function of Reynolds number. 
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Fig. 4. — Total energy spectra compensated by 
e T 2 ^ 3 fc 5 / 3 , where ex is the rate of total energy 
dissipation. Spectra for different resolution are 
compared. The dissipation cutoff wavenumber, 
kd = (ct/^ 3 ) 1 ^ 4 , is indicated by short arrows 
at the top of the plot. The Reynolds number, 
Re = u r ms/("fcf), is given in the legend. 

the physical V 2 diffusion operator and our dis- 
cretization is accurate to sixth order in space and 
third order in time. We have compared two runs 
with identical values of v and r], one with 256 3 



meshpoints and the other one with 512 3 , and 
found the same spectra (Haugen, Brandenburg, & 
Doblcr 2003). We have also compared with runs 
using double precision and found the spectra to be 
the same. The runs with up to 512 3 meshpoints 
have run for up to 80 turnover times, (Mpmsfe) -1 , 
which should be long enough to eliminate tran- 
sients. The run with 1024 3 meshpoints has only 
run for 5 turnover times, but the results are oth- 
erwise in qualitative agreement with those of the 
512 3 run. 

By comparing with nonmagnetic runs at the 
same resolution (1024 3 meshpoints) we have found 
that in hydromagnetic turbulence the bottleneck 
effect is about equally strong. This leaves us with 
the possibility that the bottleneck effect is real, 
both for hydrodynamic and for hydromagnetic 
turbulence, but that it is simply more pronounced 
in fully three-dimensional spectra, ft is there- 
fore still possible that at larger Reynolds numbers 
the true inertial range spectrum will have a k~ 5 / 3 
behavior both for kinetic and magnetic energies 
(Goldreich & Sridhar 1995). 

Alternatively, a k~ 3 > 2 spectrum might be read- 
ily explicable in terms of the Iroshnikov-Kraichnan 
phenomenology (Iroshnikov 1963; Kraichnan 1965). 
It does of course ignore local anisotropy, but more 
importantly, it predicts that the fourth order 
structure function, Si(t), scales linearly in the 
inertial range, (4 = dlnS^/ dint = 1; sec Biskamp 
(1993). To assess this possibility, we calculate the 
double-logarithmic derivative, £ p = dlnS^/dln^, 
of the unsigned structure function, S p — (|2; = ' = | p ), 
where = u ± B/^/iqP are the Elsasser vari- 
ables. The scalings of z + and z~ turn out to be 
similar, so we take the average value. 

Figure 5 shows that our data are inconsis- 
tent with Iroshnikov-Kraichnan scaling (£4 w 1.3, 
rather than 1). Instead, our data are consis- 
tent with £3 = 1 in the inertial range (Fig. 6). 
Note that linear inertial range scaling of Ss(£), 
i.e. £3 = 1, is an exact result for hydrodynamic 
turbulence (Frisch 1995). Using the extended self- 
similarity hypothesis (Benzi 1993) we plot S p (£) 
against Ss(£) to obtain a more accurate determi- 
nation of ( p for p ^ 3. Our results are compat- 
ible with She & Leveque (1994), and hence also 
with earlier forced simulations with an external 
field (Cho et al. 2002a). In particular, (2 + 1 is the 
negative slope of the total energy spectrum: we 
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Fig. 5. — Fourth order structure function for runs 
with 512 3 meshpoints. Clearly, S4 (£) is not com- 
patible with linear scaling. The inset gives the re- 
sult for 256 3 meshpoints. The scaling for transver- 
sal structure functions (dotted lines) tends to be 
better than for the longitudinal ones (solid lines). 
The statistics for the 256 3 runs is somewhat better 
than for the shorter 512 3 runs. 
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Fig. 6. — Same as Fig. 5, but for S${1). Note that 
53 {£) is compatible with linear scaling, i.e. £3 = 1 
in the inertial range. 

find (2 = 0.70, which is again in agreement with 
Kolmogorov scaling. This agrees also with sim- 
ulations of decaying helical turbulence (Biskamp 
& Miiller 2000), except that there the higher mo- 
ments have been found to be smaller. 

4. Conclusions 

For nonhelically forced hydromagnctic turbu- 
lence a resolutions of 1024 3 is necessary in order 
to begin to establish inertial range scaling. Nev- 
ertheless, the energy spectra show what we inter- 



pret as a strong bottleneck effect. This effect is 
particularly strong in the shell-integrated spectra. 
Throughout the inertial range, however, the mag- 
netic energy exceeds the kinetic energy by a factor 
of about 2.5. Both kinetic and magnetic energies 
are dominated by the spectral values at the be- 
ginning of the inertial range and independent of 
magnetic resistivity. The values (2 — 0.70 and 
(3 = 1.0 strongly favor an asymptotic fc -5 / 3 spec- 
trum. 

Our results do not support recent claims that 
the magnetic energy spectrum peaks at the resis- 
tive scale (Maron & Cowley 2001; Maron & Black- 
man 2002; Schekochihin ct al. 2002a, b). However, 
closer inspection of Fig. 3 of Maron & Cowley 
(2001) or Fig. 1 of Maron & Blackman (2002) re- 
veals that also in their cases the magnetic energy 
peaks at about 5, in agreement with our results. 

1 /2 

This is comparable to the value k{R c ~ 8 ex- 
pected based on a nonlinear closure model (Sub- 
ramanian 1999). Here, R c w 30 is the critical mag- 
netic Reynolds number for dynamo action (Hau- 
gen, Brandenburg, & Doblcr 2003), and fcf w 1.5. 
For k > 8 we find the spectrum to be Kolmogorov- 
like. 

In the interstellar medium, the magnetic Prandtl 
number is very large (Kulsrud & Anderson 1992). 
This is also the regime in which the claims regard- 
ing the peak at the resistive scale are thought to 
apply best. On the other hand, preliminary results 
suggest that even for a magnetic Prandtl number 
between 5 and 30 the magnetic energy spectrum 
still peaks at k w 5 (Haugen, Brandenburg, & 
Doblcr 2003). However, the spectrum shows now 
a possible k^ 1 tail near the viscous dissipation 
range where ohmic dissipation is still weak (Cho et 
al. 2002b; Haugen, Brandenburg, & Dobler 2003). 
In the interstellar medium, even though the mag- 
netic Prandtl number is very large, the Reynolds 
number is also very large, so the viscous cutoff 
scale is still small (~ 10 16 cm = 0.003 pc). We 
expect that in the range 0.01 . . . lOpc kinetic and 
magnetic energy spectra are parallel and show 
sca ling, but currently the numerical resolu- 
tion is still insufficient to demonstrate this. 

We thank Jason Maron, Ake Nordlund and 
Kandaswamy Subramanian for their comments on 
the paper. Use of the parallel computers in Trond- 
heim (Gridur), Odense (Horseshoe) and Leicester 
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